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OPTIMAL TIME-DEPENDENT LOWER BOUND ON DENSITY 
FOR CLASSICAL SOLUTIONS OF 1-D COMPRESSIBLE EULER 

EQUATIONS 


GENG CHEN 


Abstract. For the compressible Euler equations, even when initial data 
are uniformly away from vacuum, solutions can approach vacuum in infi¬ 
nite time. Achieving sharp lower bounds of density is crucial in the study 
of Euler equations. In this paper, for the initial value problems of isen- 
tropic and full Euler equations in one space dimension, assuming the initial 
density has positive lower bound, we prove that density functions in clas¬ 
sical solutions have positive lower bounds in the order of 0(1 + t) _1 and 
0(1 + t)~ 1 ~ s for any 0 < S <C 1, respectively, where t is time. The orders 
of these bounds are optimal or almost optimal, respectively. Furthermore, 
for classical solutions in Eulerian coordinates (y,t) £ R x [0,T), we show 
velocity u satisfies that u y (y, t) is uniformly bounded from above by a con¬ 
stant independent of T, although u y (y 1 t ) tends to negative infinity when 
gradient blowup happens, i.e. when shock forms, in finite time. 
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1. Introduction 


The compressible Euler equations in Lagrangian coordinates in one space dimen¬ 
sion are 


( 1 . 1 ) 

( 1 . 2 ) 

(1.3) 


Tt -U x = 0 

Ut+p x = 0 



where p is the density, r = p is the specific volume, p is the pressure, u is the 
velocity, e is the specific internal energy, t £ is the time and x G M is the 
spatial coordinate. The compressible Euler equations are widely used, especially in 
the gas dynamics. The classical solutions for the compressible Euler equations in 
Lagrangian and Eulerian coordinates are equivalent [lOj . 

For simplicity, in this paper, we only consider the case when the gas is ideal 
polytropic, in which 


p = K e c v t 7 with adiabatic gas constant 7 > 1, 


(1.4) 


and 
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where S is the entropy, K and c v are positive constants, c.f. [9] or [16]. For C 1 
solutions, it follows that is equivalent to the conservation of entropy [IB] : 

S t = 0 , (1.5) 

hence 

S(x,t ) = 5(x,0) = S(x). 

If the entropy is constant, the flow is isentropic, then (11.111 and m become a 
closed system, known as the p-system: 

r t -u x = 0 , 

ut + p x = 0 , 

with 

p = K r -7 , 7 > 1, (1.8) 

where, without loss of generality, we still use K to denote the constant in pressure. 

In this paper, we consider the classical solutions of initial value problems for full 
Euler equations (11.111 . (11.211 . (11.411 and (11.511 with initial data (u(x, 0), t(x, 0), S(x, 0)) 
and isentropic Euler equations (11.611 r\_/ (11.81) with initial data (u(x, 0), t(x, 0)). We 
consider the large data problem, which means that there is no restriction on the size 
of the solutions. 

Toward a large data global existence of BV solutions for the compressible Euler 
equations, which is a major open problem in the field of hyperbolic conservation 
laws, one of the main challenges is the possible degeneracy when density approaches 
zero. In fact, a solution loses its strict hyperbolicity as density approaches zero. 
See [1, .SE] for analysis and examples showing these difficulties. Therefore, the 
sharp information on the time decay of density lower bound is critical in the study 
of compressible Euler equations. Furthermore, the time-dependent lower bound on 
density for classical solutions can be used to study the shock formation and life-span 
of classical solutions. 

The study of lower bound of density for classical solutions can be traced back to 
Riemann’s pioneer paper m in 1860, in which he considered a special wave inter¬ 
action between two strong rarefaction waves. By studying Riemann’s construction, 
Lipschitz continuous examples for isentropic Euler equations ([1.611 ([LSD were pro¬ 
vided in Section 82 in [9], in which the function rnin^gR p(x, t) was proved to decay 
to zero in an order of 0(1 H-t )” 1 as t —» oo, while the initial density is uniformly away 
from zeroj^ A relative detailed discussion can be found in [5], when the adiabatic 
constant 7 = with any positive integer A f. 

Then there were many articles working on time-dependent lower bound on den¬ 
sity for general classical solutions of isentropic Euler equations (11.611 ([LSD under 
assumption that initial density is uniformly positive. For rarefactive piecewise Lip¬ 
schitz continuous solutions, for any 7 > 1, L. Lin first proved that the density has 
lower bound in the order of 0(1 +1 ) -1 in [13] by introducing a polygonal scheme. 
A breakthrough for general classical solutions happens in a recent paper a , in 
which R. Pan, S. Zhu and the author found a lower bound of density in the or¬ 
der of 0(1 + t) _4 /( 3_ '>') when 1 < 7 < 3. Using this result together with Lax’s 
decomposition in ns], Pan, Zhu and the author proved that gradient blowup of u 

4 The author thanks Helge Kristian Jenssen who first pointed out this result to him. 


( 1 . 6 ) 

(1.7) 
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and/or r happens in finite time if and only if the initial data are forward or back¬ 
ward compressive somewhere. Next, for general Lipschitz continuous solution, Pan, 
Zhu and the author in [5] improved the lower bound on density from the order of 
0(1 -|-t) -4 /( 3- T) to the optimal order 0(1 + t)~ 1 by introducing a polygonal scheme. 
The advantage of this method is that it works for not only classical solutions but 
also Lipschitz continuous solutions. And the scheme itself is of both analytical and 
numerical interest. However, the use of a polygonal scheme makes the proof very 
complex and the method seems hard to be extended to full Euler equations. Another 
result on the lower bound of density for classical solution in the order of 0(1 + f) -1 
when 7 = 3 was given by A. BressarH, where the proof relies on the study of Riccati 
equations established by Lax in m- 

For non-isentropic full Euler equations, before this paper, the only polynomial 
order upper bound of r (lower bound of p ) for general classical solution was estab¬ 
lished by Pan, Zhu and the author in [4|. More precisely, we showed density has a 
lower bound in the order of 0(1 + t) _ 4 A 3_ '>') w hen 1 < 7 < 3. 

In summary, lower bound of density in optimal order 0(1 + t ) -1 is still not 
available for isentropic Euler equations with 7 > 3 and full nonisentropic Euler 
equations with 7 > 1 , before this paper. 

In this paper, we consider classical solutions of Cauchy problems of both isentropic 
Euler equations and nonisentropic Euler equations. And we assume that initial 
density is uniformly positive. We give a short proof that density has time-dependent 
lower bound in optimal order 0(1 for isentropic Euler equations (in Theorem 

12.11) and in almost optimal order 0(1 + for any 0 < 5 < ^ for full Euler 

equations (in Theorem 12.31) in one space dimension, respectively. 

Furthermore, for classical solutions, we prove that u x (x,t) for p-system and p e u x 
for any 0 < e < ^ for full Euler equations are uniformly bounded above by a con¬ 
stant, respectively, although they are unbounded from below when gradient blowup 
happens, i.e. when shock forms. In Eulerian coordinates ( y,t ), we show for full 
Euler equations u y (y,t ) is uniformly bounded above by a constant. 

The lower bounds of density achieved in this paper can give us more precise 
estimate of life span of classical solution than those achieved in jlj and motivate us 
in searching lower bound of density for BV solutions including shock waves, which 
is a major obstacle in establishing large BV existence theory for Euler equations. 
Another interesting result on a time-dependent density lower bound for isentropic 
Euler-Poisson equations can be found in m by E. Tadmor and D. Wei. 

The rest of the paper is divided into three sections. In Section 2, we introduce 
the main results and ideas in this paper. In Section 3, we prove Theorem 12.11 for 
the p-system. In Section 4, we prove Theorem 12.31 for the full Euler equations. 

2. Main results and ideas 

Let’s first introduce some variables and notations. For Euler equations m r^j (ESD, 
we use variables 

m = e 2^17 and p = ^ 777 ^ ( 2 . 1 ) 

to take the roles of S and r. We denote the Riemann invariants 

s=u+mp and r=u—mp ( 2 - 2 ) 

2 The author knew this unpublished result through a private communication with A. Bressan. 
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respectively, and gradient variables 

a = u x + m rj x + J ^-m x r) and fi = u x - m r) x - ' 1 ^ L m x r). (2.3) 

For the isentropic Euler equations (p-system) (11.61) D, whose solutions are 
special solutions of full Euler equations (11.11) ^ (11.41) when we restrict our considera¬ 
tion on the classical solution, Riemann invariants are 


s = u + and r = u — rj 


(2.4) 


and 

a = u x + rj x = s x and /3 = u x - p x = r x . (2.5) 

The main results in this paper are listed in the following two theorems: Theorem 
IQ for p-system and Theorem 12.31 for full Euler equations. 

Theorem 2.1. Let (u(x, f), r(x, t)) be a C 1 solution of the isentropic Euler equa¬ 
tions (EH) \1.8\) in the region (x,t) € M x [0,T), where T can be any finite positive 
constant or infinity. Assume ri(x,0), r(x,0) > 0, p(x, 0) = l/r(x,0), a(x,0) and 
(3(x, 0) are all uniformly bounded, where a and j3 take the form in (12.5p . 

Let M be an upper bound of a{x, 0) and f3(x, 0), i.e. 



max < a 
xGM l 

(x,0),/3(x,0)} < M 

(2.6) 

then 

max 

^a(x,t), f3(x,t)^ < M . 

(2.7) 


(:c,t)EMx [0,T) 

This gives 

max {rt 

(#,£)E Rx [0 ,T) 

) = max | < M 

’ (x,f)£Rx [0,T) k J 

(2.8) 


by 12.51) and fZ2P. Hence, there exist positive constants M\ and M 2 independent of 
T such that 

M\ 

mm p(x,t) > ■ 2.9) 

x M 2 + t 

The key step in the proof of Theorem 12.II is to prove (|2.7I) . In fact, suppose (I2.7jl 
is correct, then by the conservation of mass (11.61) and (12.5|) . we can easily prove 

AMD : 

r t = u x = ^(a +{}) < M (2.10) 

which directly gives (12.91) . together with the initial condition. To prove (12.71) . we 
need to study the characteristic decomposition established by Lax in Em The key 
idea is to find an invariant domain on a and (3. 

One conclusion that we can draw from (12.61) ^ (12.71) is that although the variables 
a and fi might increase along forward and backward characteristics, respectively, the 
function max ie R{a(r, t), (3{x, f)} is not increasing with respect to t, which means 
that the maximum rarefaction of classical solution is not increasing. This result can 
be easily seen from the fact that (12.71) is still correct if we change 0 in (12.61) into any 
t* E (0,f). 
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Remark 2.2. Under assumptions in Theorem \2.1l in Eulerian coordinates ( y,t ), 
the inequality 12.81) gives that smooth solutions in the region (y, t) E Kx [0, T) satisfy 

max { — \<M, (2.11) 

(y,t)mx[0,T) 1 P J 

where M is the constant given in m, because pu x (x,t ) = u y (y,t). See [XB| for 
the transformation between Eulerian and Lagrangian coordinates. 

Since p is uniformly bounded above, which can be easily proved by the fact that 
Riemann invariants s and r are initially bounded and are constant along forward 
and backward characteristics, respectively, we know 

max \u v \ < M , 

(y,t)eRx[0,T) 1 J 

for some constant Ad independent of T. 


Then we consider the full Euler equations. 


Theorem 2.3. Let (u(x, t), r(x, t), S(x)) be a C l solution of full Euler equations 
mn> r*j (HI in the region (x,t) € I x [0, T). Here, T can be any finite positive con¬ 
stant or infinity. Assume that initial data u(x, 0), r(x,0) > 0, p(x, 0) = l/r(x,0) ; 
S(x), S'(x), a(x,0) and f3(x, 0) are all uniformly bounded and total variation of 
S(x) is finite, where a and /3 satisfy A 2. 31) . Then, for any 


there exists constant Nq independent of T, such that 


max 

(#,t)EMx [0 ,T) 



max 

(a;,t)GRx [0,T) 




( 2 . 12 ) 


and there exist positive constants N\ and N 2 independent of T such that 

• (213 > 

where 5 = > 0 . 


We first prove a result in Lemma 14.41 taking the similar role as <\2.7\) in Theorem 
12.91 In fact, we find uniform bounds on gradient variables p £ a and p £ (3, using which 
we can easily prove (12.12(1 by (|2.31) and (jl.lf) : 

p £ r t = p £ u x = ^{p £ a + p £ /3) < Constant, 

then show (12.1311 . The reason why we use p £ a and p £ /3 instead of a and /3 is to 
control the lower order terms in the Riccati equations produced by the varying 
entropy. The proof of Theorem 12.31 also relies on the uniform constant upper bound 
of density established in [ 8 ] by R. Young, Q. Zhang and the author for classical 
solutions when total variation of initial entropy is finite. 


Remark 2.4. Under assumptions in Theorem 1 2.31 in Eulerian coordinates ( y,t ), 
the inequality i2.8f) gives that the classical solution in the region ( y , t) € R x [0, T) 
satisfies 


max 

(j/,t)eMx[0,T) 



<N 0 . 
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Since p is uniformly bounded above under assumptions in Theorem ) 2. 31 we know 


max 

(j,i)£Rx [0,T) 


jtiyj < N 0 , 


for some constant No independent of T. 

See [16] for the transformation between Eulerian and Lagrangian coordinates. 
Since this result is a local result, we only need to assume that initial entropy is 
locally BV. 


One direct application of Theorem 12.31 is that one can use (12.131) to improve the 
life-span estimate established in [4J when 1 < 7 < 3 which depends on the time- 
dependent lower bound of density. We leave this to the reader. 


3. Lower bound of density for p-system: The proof of Theorem 12.11 


We first introduce the characteristic decompositions for C 1 solution of p-system. 
For any classical solution for (11.61) ^ (11.81) . the Riemann invariants s and r in (12.41) 
are constant along forward and backward characteristics, respectively, 


d + s = 0 and <9_r = 0 


(3.1) 


with 


d + = dt + cd x and d- = dt — cd x 


and wave speed 

c = V~Pr = y/K 71 - 2 “ . 

Furthermore, gradient variables a = s x and fd = r x defined in (12.51) satisfy the 
following Riccati equations. 


Proposition 3.1. [2] The classical solution in m 11.81) satisfy 


d + a = ki{a/3 — a 2 } , (3-2) 

and 

d-/3 = ki{a/3 — /3 2 } , (3.3) 


where 


k\ = 


(7+1 )K 
2 ( 7 - 1 ) 



(3.4) 


where K c is a positive constant given in m ■ The function p > 0 is defined in 

m- 

Equations (13.21) and (13.31) are special examples of Lax’s decompositions in [12] for 
general hyperbolic systems with two unknowns. See detailed derivation of (13.21) and 

m> in El- 


Remark 3.2. The idea for the proof of (|2.7D can be seen from FigureUi 


Before the proof, we remark on one fact that p , p, c and k\ are all bounded above 
by some constants if assumptions in Theorem 12.11 are satisfied. This can be easily 
obtained by (13.11) , which says that s and r are constant along forward and backward 
characteristics. As a consequence, p , p, c and function k\ are all uniformly bounded 
from above. Denote 

K\ = max ki(x,t ), (3.5) 

(z,t)6Kx[0,T) 

where K\ is a constant only depending on 7 and initial condition. 
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Figure 1. max{a:,/3} < M is an invariant domain. Note: a (or j3) 
might increase along characteristic. 

Proof of Theorem 12.11 We first prove (12.71) by contradiction. Without loss of 
generality, assume that a(xo,to) = M at some point (xo,to). See Figure El 



Figure 2. Proof of Theorems 12.II and 12.31 

Because wave speed c is bounded above, we can find the characteristic triangle 
with vertex (xo,to) and lower boundary on the initial line t = (1. denoted by 12. 

Then we can find the first time t± such that a = M or p = M in 12. More 
precisely, 

max (a(x,t), f3(x,t)) < M, 

and a(xi,t\) = M or/and /3(x\,ti) = M for some (xi,ti) € Q. Without loss of 
generality, still assume a(xi,t\) = M. The proof for another case is entirely same. 
Let’s denote the characteristic triangle with vertex as fii G 12, then 

max (a(x,t), (3(x,t)) < M, (3.6) 

(x,*)€«!,*<*! V V 

and cc(xi,ti) = M. By the continuity of a, we could find a time 1 2 G [0, t±) such 
that, 

a(x,t) > 0, for any (x,t) G 12i and t > ■ (3-7) 

Next we derive a contradiction. By (13.21) . (13.51) and (|3.6j) ~ (j3.7j) . along the forward 
characteristic segment through (xi,fi) when t 2 < t < t\, 

d + a = ki{a/3 — a 2 } < Ki{Ma — a 2 } 
which gives, through integration along characteristic, 

(M^ < *1 dt 

^ W ln M—a(t) — if M-a(t 2 ) + _ • 
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As t —> 1 1 —, left hand side approaches infinity while right hand side approaches a 
finite number, which gives a contradiction. Hence we prove that (12.71) is correct, i.e. 
a and (3 are uniformly bounded above. Then by the conservation of mass m and 
(12.51) . we have (12.101) then (12.81) . which directly gives (12.91) . Hence we complete the 
proof of Theorem 12.11 

Remark 3.3. Theorem \2.1\ can he extended to the case with general pressure law 
P = p{t) with p T < 0, p TT > 0 and some other suitable conditions on p. We 
leave this to the reader and refer the reader to [ 6 j for the Riccati equations and the 
definitions of a and (3. For full Euler equations, the extension of Theorem 1 2. .VI to 
general pressure law is still not available because the current result on uniform upper 
bound of density is only available for 7 -law pressure. 


4. Full compressible Euler equations 


4.1. Equations and coordinates. We first introduce some notations and existing 


equations for C 1 solutions of full Euler equations (11.11) m- Recall we use new 
variables m and 7 to take the roles of S and r, respectively: 


s 

m = e 2c >' 


(4.1) 


and 



(4.2) 


Without confusion, we still use c to denote the nonlinear Lagrangian wave speed for 
full Euler equations, where 


7+1 £ 

2 g2 c v . 


V~Pt = V k t t 


(4.3) 


c = 


The forward and backward characteristics are described by 


and we denote the corresponding directional derivatives along these characteristics 


by 



respectively. 

It follows that 


2 


r = K t 7- 1 , 

2 - 23 - 

p = K p m 77 T - 1 , 

7+1 

c = c(? 7 , m) = K c mr) t - 1 . 


(4.5) 


with positive constants 



K V -.= KK~^, and K c := sfK^ f K T 2 , (4.6) 


K P = ^K c and K T K C = ^ . 


so that also 


(4.7) 
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In these coordinates, for C 1 solutions, equations (11.11) — (I1.4D are equivalent to 

rit + — u x = 0, (4.8) 

m 

P 

ut + mcrjx + 2— m x = 0 , (4.9) 

m 

m t = 0 , (4-10) 

where the last equation comes from (11.51) . which is equivalent to (11.31) . c.f. fl6j . 


Note that, while the solution remains C 1 , m = m{x) is given by the initial data and 
can be regarded as a stationary quantity. 

Recall that we denote the Riemann invariants by 

r:=u — mr] and s:=u + mrj. (4-11) 

Different from the isentropic case (m constant), for general non-isentropic flow, s 
and r vary along characteristics. Also recall we denote gradient variables 

a = u x + mr) x + (4.12) 

P = u x - mrj x - Tilm x r) , (4.13) 

which satisfy the following Riccati equations. See detailed derivation in |2j. 
Proposition 4.1. |2] The classical solutions for El El satisfy 

d + a = fci{l’2(3a + /?) + a/3 — a 2 }, (4-14) 

and 

d-P = k\{— ^ 2(0 + 3/3) + a(3 — /3 2 }, (4.15) 

where 

fci = k 2 = ^^V m x- (4-16) 

Proposition 13.11 is in fact a corollary of Proposition 14.11 for the isentropic case in 
which m x = 0. 

4.2. Uniform upper bound on density. In this part, we review a result on the 
uniform upper bounds of |ri| and p established by the author, R. Young and Q. 
Zhang in [7], for later references. 

In this section, we always assume all initial conditions in Theorem l2.3l are satisfied. 
So that 

1 f + °° r + °° I m'(x)\ 

V:= - / \S'(x)\ dx = / dx < oo, (4.17) 

2 c v J.oo J-cc m{x) 

while also, by (ED, 

0 < Ml < m(-) < Mu , (4-18) 

for some constants Ml and Mu. Also there exist positive constants M s and M r , 
such that, in the initial data, 

|s(-, 0)| < M s and |r(-,0)| < M r . (4-19) 

In the following proposition established in |7J, |w| and p are shown to be uniformly 
bounded above. 
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Proposition 4.2. [7j Assume all initial conditions in Theorem \2.A are satisfied. 
And assume system El El has a C 1 solution when t € [0, T), then one has the 
uniform bounds 

\u(x,t)\ < Mu *7 and rj(x,t ) < ^ 1 + Ml t ' " 1 , (4.20) 

where 

Li := M s + V M r + V (V M s + V 2 M r ) e ^, 

L 2 := M r + VM S + V(V M r + V 2 M s ) +', 

and 



Constants L\ and L 2 both clearly depend only on the initial data and 7 . Here T can 
be any positive number or infinity. And the bounds are independent of T. 

4.3. Proof of Theorem 12.31 Similar as Theorem 12.II for p-system, the key idea is 
still to get the uniform upper bound of some gradient variables measuring rarefac¬ 
tion. 

However, we cannot directly get the uniform upper bound of a and /3. In fact, 
comparing to (13.21) ^ (13.31) . equations (I4.14IW(|4.15I) include some first order terms in 
the right hand side. In order to cope with them, we introduce some new gradient 
variables 

2e 2e 

a £ = 7 T - 1 a and (3 e = rj^- 1 /3. (4-21) 

Using (14.81) . we have 

d+rj = rj t + aix = ~—u x + cq x = - 1 —K c r]^- 1 m x , 

m 1 

and 

C 7+1 _H 27 

<9-7 = rj t — cq x = — -u x - crj x = - K c ryy~ 1 a + T—Kciyi- 1 m x , 
m 7 

then it is easy to prove the next lemma by Proposition 14.11 


Lemma 4.3. The classical solutions in El El satisfy 



d+ot £ = ki £ {M3a £ - 4 ea E + /3 £ ) + (1 - ^-)a £ P £ - a%j , 

(4.22) 

and 

d-/3 £ = k l£ {-k 2£ (a £ + 3f3 £ - 4ef3 £ ) + (1 - -)a £ /3 £ - $}, 

(4.23) 

where 

fei — ^ +1 )K C T (l-s) 7 _ 7-1 l+Tzje 

~ 2(7-1) " 1 7(7+1) '/ rn xi 

(4.24) 

and 

1 

0 < £ < -. 

4 

(4.25) 

Note, 

for any C 1 solutions in (x,t) £ lx [ 0 , T) satisfying initial conditions in 


Theorem 12.31 using Proposition 14.21 for any e satisfying (14.251) . we know \ki £ (x,t)\ 
and \k 2e {x,t)\ are both uniformly bounded above: 


\ki e (x,t)\ < Ki and \k 2e (x, t)\ < K 2 , 


(4.26) 
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where constants K\ and K 2 only depend on initial conditions and 7 but independent 
of e. 

Next we give the key lemma which will be proved later. 


Lemma 4.4. Suppose the initial conditions in Theorem \ 2. S\ are satisfied. For any 
e satisfying (14.251) . let N be an upper bound of a e (x, 0) and f3 e (x,0), i.e. 

max ia £ (x, 0), fi e (x, 0)1 < N (4.27) 


where constant N also satisfies 

jV> m ax{! 

7+1 

then 

max \ a E (x,t), BJx,t) \ < N. 
(+,i)eRx[0,T) l v ’ v 

mB- 


(4.28) 

(4.29) 


Proof of Theorem 12.31 We only have to show Lemma 14.41 In fact, if Lemma 14.41 
is proved, then by the conservation of mass (jl.ll) and definitions of a £ and fi £ in 
(|4.21j) and (14.121) ^( 14. 131) . we have 

2e 2e 1 

7 7_ l T t = T)"!- 1 u x = -{a £ + fi £ ) < N 

which gives that, by (14.21) . t = 1/p and initial density has positive lower bound, 
there exists positive constants N\ and A7 2 , such that 


P > ( 


iVi 


n 2 + t 


1+5 


where 


Then it is easy to see that all results in Theorem 12.31 are correct. 

Now we prove Lemma l4.4l bv contradiction. We still use Figure [2j Without loss 
of generality, assume that a £ (xo,to) = N, at some point (xo,to)- 

Because wave speed c is bounded above, we can find the characteristic triangle 
with vertex (xo,to) an d lower boundary on the initial line t = 0 , denoted by Q. 

Then we can find the first time t\ such that o. £ = N or fj £ = N in O. More 
precisely, 

max \a £ (x,t),fi e (x,t)) < N, 

and a e (xi,ti) = N or/and f3 £ (xi,ti) = N for some (xi,t\) € fl. Without loss of 
generality, still assume a e (xi,t\) = N. The proof for another case is entirely same. 
Let’s denote the characteristic triangle with vertex (xi,ti) as fii G Cl, then 


max 

(. x,t)eQi,t<ti 


a £ (x,t),/3 £ (x,t)^ 


<N, 


and a £ (xi,t\) = N. 

Then we divide the problem into two cases: 


I. N > (3 e {xi,ti) > ; 

II. P e (xi,tl) < -f- 
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In case I, by the continuity of a £ and ft £ and our construction, we can find a time 
t ,2 E [0, ti) such that, 

N 

— < a E (x,t ) < N and |/3 e | < N , for any (x,t) € Hi and t 2 < t < t\. 

(4.30) 

Then using (14.221) , (14.261) , (|4.28[) and (14.301) , along the forward characteristic segment 
through (xi,fi), when t 2 < t < t\, we have 


d + a £ < k le (l - (a £ /3 £ - a^) < (JVa £ - a^) 


with 


K\ = Kfl (1 - 


4e 


y + 1 

which gives, through integration along characteristic, 

dae < K\dt 


( N-a e )a e 


1 


In 


a e(t) 


< 


1 


ln N-tLs + R i(t- h) 


N N—a e (t) — N N—a e (t 2 ) 

As t —> t\— , left hand side approaches infinity while right hand side approaches a 
finite number, which gives a contradiction. 

In case II, by the continuity of a e , we could find a time ts E [0, ti) such that, 

N N 

— < a £ (x,t) < N and (3 e (x,t) < -, for any (x,t) E flj and t% < t < t± . 

(4-31) 

which gives, by (14.281) . 

{k 2e + (1 - )a e ) f3 e < 0. 

Hence by (14.281) . (14.251) and (14.311) . we have 

d + a £ < ki e {/C2e(3 - 4e)a £ - a £ } < 0 . 


As a consequence, a £ decreases on t along the forward characteristic line through 
(xi,ti), when t% < t < ti, which contradicts to that a £ (xi,t\) = N while a £ (x,t ) < 
N when (x,t) E Hi and < t < t\. Hence Lemma 14.41 is proved. This completes 
the proof of Theorem 12.31 
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